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Abstract 

In  this  paper,  we  describe  direct-sequence  code-division  multiple-access  (DS-CDMA)  systems  with  quadriphase-shift  keying 
in  which  channel  estimation,  coherent  demodulation,  and  decoding  are  iteratively  performed  without  the  use  of  any  training  or 
pilot  symbols.  An  expectation-maximization  channel-estimation  algorithm  for  the  fading  amplitude,  phase,  and  the  interference 
power  spectral  density  (PSD)  due  to  the  combined  interference  and  thermal  noise  is  proposed  for  DS-CDMA  systems  with 
irregular  repeat-accumulate  codes.  After  initial  estimates  of  the  fading  amplitude,  phase,  and  interference  PSD  are  obtained  from 
the  received  symbols,  subsequent  values  of  these  parameters  are  iteratively  updated  by  using  the  soft  feedback  from  the  channel 
decoder.  The  updated  estimates  are  combined  with  the  received  symbols  and  iteratively  passed  to  the  decoder.  The  elimination  of 
pilot  symbols  simplifies  the  system  design  and  allows  either  an  enhanced  information  throughput,  an  improved  bit  error  rate,  or 
greater  spectral  efficiency.  The  interference-PSD  estimation  enables  DS-CDMA  systems  to  significantly  suppress  interference. 
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I.  INTRODUCTION 

In  mobile  communication  systems,  the  wireless  channel  induces  random  amplitude  and  phase  variations  in  the  received 
data,  with  the  possible  addition  of  time-varying  interference  from  co-channel  users.  For  this  reason,  the  accuracy  of  channel 
state  information  (CSI)  at  the  receiver  is  critical  for  coherent  detection  and  demodulation.  A  number  of  methods  have  been 
proposed  for  estimation  of  CSI,  all  of  which  fall  within  the  broad  categories  of  either  pilot-assisted  or  blind  algorithms.  Current 
and  next-generation  cellular  protocols  such  as  W-CDMA  (Wideband  Code  Division  Multiple  Access)  and  3GPP  LTE  (Third 
Generation  Partnership  Project  Long-Term  Evolution)  specify  the  use  of  pilot-assisted  channel  estimation  (PACE)  C D-  Pilot 
symbols  or  training  sequences  are  known  symbols  either  multiplexed  with  or  superimposed  onto  the  transmitted  data  in  the  time 
or  frequency  domain,  with  the  associated  disadvantage  of  a  loss  in  spectral  and/or  power  efficiency.  Moreover,  superimposed 
PACE  is  degraded  at  low  signal-to-noise  ratios,  and  multiplexed  PACE  is  unsuitable  for  fast-fading  channels  with  a  coherence 
time  shorter  than  the  pilot-symbol  transmission  rate  El,  E). 

Blind  channel-estimation  methods  offer  an  alternative  approach  that  avoids  the  implementation  cost  of  pilot  symbols  J4|. 
Blind  methods  typically  use  second-order  statistics  of  the  received  symbols  for  CSI  estimation,  with  shortcomings  such  as 
increased  complexity,  slow  convergence  times,  and  channel-phase  ambiguity  0-  In  addition,  the  received  interference  power 
spectral  density  (PSD),  which  is  due  to  both  the  thermal  noise  and  the  time-varying  interference,  is  usually  not  estimated  in 
the  literature  spanning  both  PACE  and  blind  CSI  estimation.  The  accuracy  of  the  interference-PSD  estimation  is  known  to  have 
a  significant  impact  on  turbo-principle  (iterative)  detection  techniques  as  well  as  turbo  and  low-density  parity-check  (LDPC) 
channel  decoding  |[6l,  171. 

The  expectation-maximization  (EM)  algorithm  offers  a  low-complexity  iterative  approach  to  optimal  maximum-likelihood 
detection  and  estimation  0,  El-  A  substantial  body  of  literature  can  be  found  on  EM-based  techniques  for  data  detection, 
multiuser  detection,  channel  estimation,  or  a  combination  of  the  latter.  A  few  representative  examples  are  listed  next.  A  recursive 
estimation  of  the  fading  channel  amplitude  was  proposed  in  flOl.  Iterative  receivers  with  EM-based  fading-amplitude  and  data 
estimation  using  pilot  symbols  for  LDPC-based  space-time  coding  and  space-time  block-coded  orthogonal  frequency-division 
multiplexing  (OFDM)  were  studied  in  1111  and  lfl2l.  respectively.  Joint  multiuser  detection  and  channel/data  estimation  for 
uplink  code-division  multiple  access  (CDMA)  was  studied  in  lfT3l-lfl6l.  In  [T7|],  iterative  EM  estimation  and  turbo  coding 
were  studied  assuming  noncoherent  frequency-shift  keying  modulation  and  demodulation,  which  is  well-known  to  be  less 
power-efficient  than  coherent  modulation  QjO. 
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In  Oil,  an  EM  estimation  approach  for  turbo-coded  single-user  iterative  CDMA  receivers  with  binary  phase-shift  keying 
was  considered.  In  [20)  and  [21],  the  authors  replaced  turbo  codes  with  regular  LDPC  codes;  however,  II 1 911-112 11  all  featured 
as  much  as  a  9.1%  pilot-symbol  overhead  for  channel-amplitude  and  interference-PSD  estimation.  Recently,  EM-based  channel 
and  noise  estimation  techniques  were  proposed  in  l22l  and  f23l  for  multiple-antenna  systems  with  convolutional  coding  and 
as  much  as  a  10%  pilot-symbol  overhead  for  initial  channel  estimation. 

Although  the  primary  role  of  pilot  symbols  in  most  cellular  standards  is  channel  estimation,  pilot  symbols  often  play  a 
secondary  role  in  cell,  frame,  or  symbol  synchronization.  However,  alternative  methods  of  synchronization  may  be  used  when 
pilot  symbols  are  unavailable  1 18),  (24),  [25] .  In  this  paper,  a  doubly  iterative  direct-sequence  CDMA  (DS-CDMA)  receiver 
featuring  iterative  EM  channel  estimation  and  iterative  detection  and  decoding  without  any  pilot  symbols  is  presented.  The 
general  form  of  the  proposed  blind  channel  estimator  provides  fading-amplitude,  phase,  and  interference-PSD  estimates  in  both 
single-user  and  multiuser  environments,  therefore  offering  an  alternative  to  the  methods  proposed  in  j26l  and  )27l  to  rectify 
the  phase  ambiguity  of  blind  channel  estimate^  The  special  case  of  EM  channel  estimation  with  perfect  phase  information 
at  the  receiver  (e.g.,  by  means  of  a  phase-locked  loop)  is  also  considered.  The  proposed  iterative  receiver  is  capable  of  using 
higher-order  modulations  such  as  M-PSK  and  M-ary  quadrature  amplitude  modulation  (M-QAM),  although  quadriphase-shift 
keying  (QPSK)  is  demonstrated  in  this  work  for  simplicity.  In  addition,  the  proposed  system  uses  irregular  repeat-accumulate 
(IRA)  codes  instead  of  regular  LDPC  codes  for  lower  complexity  m-m. 

The  paper  is  organized  as  follows.  Section  II  describes  the  system  transmitter  and  receiver  models  including  coding, 
modulation,  and  spreading,  as  well  as  fading-channel  parameters.  Section  III  summarizes  the  proposed  EM-based  estimation 
process  that  uses  soft  feedback  from  the  channel  decoder.  Section  IV  presents  the  proposed  blind  method  for  the  initial  CSI 
estimation  and  the  possible  trade-offs  vis-a-vis  PACE.  Section  V  shows  simulation  results,  and  Section  VI  offers  conclusions. 

A  word  on  notation:  lowercase  boldface  is  used  to  represent  vectors,  while  uppercase  boldface  represent  matrices.  E  denotes 
the  statistical  expectation,  (-)T  is  the  matrix  transpose,  *  is  the  complex  conjugate,  and  [x\  is  the  largest  integer  smaller  than 
x. 


II.  SYSTEM  MODEL 


Additional  block  required  for  PACE  system 


Fig.  1 .  DS-CDMA  transmitter  with  QPSK  modulation. 


Fig.  1  shows  the  block  diagram  of  a  dual  quaternary  DS-CDMA  transmitter  ||3D  consisting  of  a  channel  encoder,  QPSK 
modulator,  and  a  direct-sequence  spreading  generator  that  multiplies  orthogonal  chip  sequences  p  /,■  and  p/  with  the  in-phase 
and  quadrature  modulator  inputs.  The  input  to  the  encoder  in  Fig.  1  is  a  binary,  independent,  identically  distributed  data  block 
of  length  K,  which  is  denoted  by  m  =  [m(l), . . .  ,m{K)\,  m(ibit)  £  [1,0]. 

'In  (26],  two  different  PSK  modulations  are  used  on  adjacent  OFDM  subcarriers  to  resolve  the  phase  ambiguity  under  slow  frequency-selective  fading.  A 
short  pilot  sequence  is  used  in  E3  to  recover  the  channel  phase,  making  it  semi-blind  in  nature.  More  importantly,  the  interference-plus-noise  PSD  is  not 
estimated  in  123  and  ED. 
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A.  Encoding,  Modulation,  and  Spreading 

Each  1  x  K  message  vector  m  is  encoded  into  a  1  x  TV  codeword  b  =  [6(1), . . . ,  b(N)]  using  a  systematic,  extended  IRA 
code  |28|.  IRA  codes  offer  a  combination  of  the  linear  complexity  of  turbo  encoding  and  the  lower  complexity  of  LDPC 
decoding  without  compromising  on  performance. 

The  (TV,  K)  IRA  code  is  constructed  following  the  methodology  proposed  in  12911,  where  the  IRA  code  parameters  were 
designed  for  use  on  a  burst-erasure  channel  with  additive  noise,  which  was  shown  to  be  a  good  surrogate  for  Rayleigh  fading 
channels.  IRA  codes  can  be  considered  to  be  a  subset  of  low-density  parity-check  codes  and  therefore  may  be  represented  by 
a  Tanner  graph  l30l.  Let  A  (a:)  =  A iX1-1  and  p  (x)  =  Yl'!'  represent  the  variable-node  and  check-node  degree 

distributions  of  the  code’s  Tanner  graph,  with  (dv,dc)  being  the  maximum  variable  and  check  node  degrees,  respectively.  Using 
density  evolution,  for  (dv  =  8  ,dc  =  7)  we  obtain  the  following  good  choices  [129] : 


A(x) 

=  0.00008  +  0.31522a;  +  0.34085a;2  +  0.0.06126a;6 
+  0.28258a;7 

p{x) 

=  0.62302a;5  +  0.37698a;6. 

(1) 

The  (TV  —  K)  x  TV  IRA  parity-check  matrix  can  be  represented  as  H  =  [Hi  |  H2],  where  sub-matrix  H2  is  a  (  TV  —  K)  x 
(TV  —  K)  dual-diagonal  matrix,  and  Hi  is  a  randomly-generated  (TV  —  K)  x  K  sparse  matrix  constmcted  such  that  H  has  the 
degree  profile  of  ([!]).  The  K  x  TV  systematic  generator  matrix  G  is  then  given  by  G  =  [I k  |  H|H2  t]  . 

For  the  simulations  in  Section  V,  Gray-labeled  QPSK  is  used  with  2  encoded  bits  mapped  into  a  modulation  symbol 
x{k)  £  { ±  1 .  ±j } ,  k  =  1 .....  -A.  Although  QPSK  is  assumed,  the  analysis  and  simulation  is  easily  extended  to  M-QAM. 
Parallel  streams  of  code  bits  are  each  spread  using  a  Gold  sequence  with  spreading  factor  g  chips/code  bit  before  rectangular 
pulse-shaping  that  produces  the  real  and  imaginary  components  of  x(k),  i.e.,  Xn(k)  =  Re  (x  (k))  and  a :j  (k)  =  Im  (a;  (fc)).  In 
practice,  an  intermediate  frequency  is  used  before  the  carrier  frequency  upconversion,  but  the  upconversion  from  baseband  to 
the  intermediate  frequency  is  omitted  for  clarity  in  Fig.  1. 

No  channel  interleaving  is  applied  to  the  IRA  code  due  to  the  inherent  interleaving  characteristics  of  the  IRA  code  itself. 
This  is  because  the  IRA  code  can  be  alternatively  represented  as  a  repetition  code  concatenated  with  a  convolutional  encoder 
(accumulator)  with  an  interleaver  between  them.  The  interleaver  is  embedded  within  the  sub-matrix  Hi  in  the  Tanner  graph 
representation  of  IRA  codes. 


B.  Channel  Model 


For  multiple-access  interference  (MAI)  environments,  the  channel  coefficients  are  generated  using  the  Jakes  correlated  fading 
model.  The  flat-fading  assumption  is  valid  when  the  information  bit-rate  is  low,  e.g.,  100  kb/s  as  usually  considered  in  this 
paper,  since  the  multipath  delay  spread  in  a  typical  cellular  environment  is  about  10  ps,  which  is  negligible  compared  to  the 
symbol  duration.  For  completeness,  the  proposed  system  and  analysis  are  extended  to  include  frequency-selective  channels  by 
including  multipath  components  with  delays  exceeding  a  chip  duration  and  using  Rake  receivers  fl8l.  BTI.  as  described  in 
Section  IV-EI  Each  codeword  or  frame  of  TV  code  bits  is  divided  into  two  different  types  of  subframes  or  blocks.  One  block 
size  is  set  equal  to  the  npB  code  bits  over  which  the  fading  amplitude  is  assumed  to  be  constant.  The  other  block  size  is  set 
equal  to  tiib  code  bits  over  which  the  interference  level  is  assumed  to  be  constant. 

Each  frame  comprises  TV/2  QPSK  code  symbols  and  Ng/2  spreading-sequence  chips  for  each  QPSK  component.  The  fading 
coefficient  associated  with  spreading-sequence  chip  c  of  either  p  /,■  or  p /  is 


C 


|_c/(nFB3)J 


c  = 


Ng 

2 


(2) 


where  Es  is  the  average  energy  per  QPSK  symbol,  a  is  the  fading  amplitude  with  E  [a2]  =  1,  and  </>  is  the  unknown 
fading-induced  channel  phase. 


C.  Iterative  Receiver  Structure 

Fig.  2  shows  a  block  diagram  of  the  proposed  dual  quaternary  iterative  receiver.  The  received  signal  is  downconverted, 
passed  through  chip-matched  filters,  and  despread  by  a  synchronized  spreading-sequence  generator  in  each  branch,  with  the 
downconverter  and  synchronization  devices  [  1_8;|  omitted  in  Fig.  2  for  clarity.  Self-interference  between  the  spreading  sequences 
of  the  desired  user  is  negligible  because  accurate  synchronization  is  assumed  at  the  receiver.  Let  TVo/2  denote  the  two-sided 
PSD  of  the  Gaussian  noise.  For  the  flat-fading  scenario,  the  complex  envelope  of  the  desired  user  at  the  kth  symbol  time  with 
active  MAI  can  be  written  as 

lk/npBi  x(k)  +  nint{k)  +  n(k),  1  <  k  <  - 


y{k)  =  c 


(3) 
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Additional  block  required  for  PACE  system 


yi(k) 
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Fig.  2.  Iterative  DS-CDMA  receiver. 


where  x(k)  is  the  complex  transmitted  code  symbol  of  the  desired  user,  n(k)  is  a  complex  zero-mean  circularly  symmetric 


Gaussian  noise  sample  with  E 

EQ. 


|rifc|  =  N0,  and  nmt(k)  is  the  interference  at  the  demodulator  due  to  interfering  users  lfl8l. 


The  time-varying  MAI  is  assumed  to  be  generated  by  interfering  users  with  a  structure  identical  to  the  desired  user,  albeit 
the  spreading  sequences  differ  and  the  fading  coefficients  are  independent.  The  despreading  in  the  receiver  tends  to  whiten  the 
interference  PSD  over  the  code-symbol  passband,  and  the  subsequent  filtering  tends  to  produce  a  residual  interference  with  a 
Gaussian  distribution.  Thus,  the  interference  PSD  due  to  the  combined  interference  and  thermal  noise  is  modeled  as  additive 
Gaussian  noise  with  a  two-sided  PSD  Iq/2  that  is  constant  over  each  block  of  uib  code  bits  but  varies  from  block-to-block. 
This  model  enables  the  derivation  of  an  EM  estimator  for  Iq  that  is  used  in  the  demodulator  metric  and  leads  to  the  suppression 
of  the  interference. 

A  receiver  iteration  is  defined  as  a  fixed  number  of  decoder  iterations  followed  by  internal  EM  iterations  in  the  channel 
estimator  of  Fig.  2,  and  then  a  single  demodulator  metric  generation.  Let 


i  denote  the  index  for  the  internal  EM  iteration  ,  i  =  1 .... ,  imax\ 

j  denote  the  index  for  the  closed-loop  receiver  iteration,  j  =  1, . . .  ,jmax- 

Let  0$  =  represent  the  estimates  of  the  fading-coefficient  and  interference-PSD  parameters  at  the  ith  EM 

iteration  during  the  jth  overall  receiver  iteration.  EM  iterations  commence  after  the  initial  channel  estimation  and  decoding, 
which  is  obtained  while  the  switch  in  Fig.  2  is  set  to  position  1.  The  subsequent  receiver  iterations  are  performed  while  the 
switch  is  set  to  position  2  in  order  to  refine  the  initial  channel  estimate  with  the  aid  of  soft  feedback  sf  \  fi  =  1,  2, 3, 4  from 
the  channel  decoder. 


III.  EM  ALGORITHM 

Theoretically,  the  maximum-likelihood  CSI  estimator  0  can  be  obtained  from  a  received  data  vector  y  =  [y(l ), . . . ,  y(Ni)] 
of  N\  code  symbols,  referred  to  as  the  incomplete  data,  by  maximizing  the  conditional  log-likelihood  function: 

0  =  argmaxln/(y  |  0).  (4) 

0 

However,  the  computation  of  this  equation  is  virtually  prohibitive  in  practice  since  its  complexity  increases  exponentially  with 
the  observation  window  size.  In  the  EM  algorithm,  the  expectation  of  the  conditional  log-likelihood  of  the  complete  data 
z  =  (x,  y)  is  iteratively  maximized  with  respect  to  0,  where  expectation  is  taken  with  respect  to  x  given  y  and  a  previous 
estimate  of  0. 

The  conditional  probability  density  function  (pdf)  of  z  can  be  written  as 

/( z  |  9)  =  /(X,  y  |  0)  =  f(y  |  X,  0)/(x  I  0)  =  /( y  I  X,  0)/(x)  (5) 

where  the  last  equality  is  from  the  independence  of  the  transmitted  signal  vector  x  and  the  CSI  parameter  0.  Thus, 

In  /  (z  |  0)  =  In  /  (y  |  x,  0)  +  In  /  (x) . 


(6) 
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Since  the  symbols  are  independent  and  circularly  symmetric  Gaussian  noise  and  interference  are  assumed,  the  conditional  pdf 
/  (y  |  x,  0)  is 

1 


/( y  |x,0)  = 


(  {\y{k)  —  Cx{k)\2 

exP  -  - 7 - 


(7) 


k= 1 


*0 


Therefore,  as  |x(fc)| 2  =  1  Vfc, 


JVi 


ln/(y|x,6/)  =  -iVi -ln(/o)  -  —  ^][|y(fc)|2  +  ICI5 


fc=i 


-2Re(y*(fe)C'*(fc))] 


(8) 


where  an  irrelevant  constant  has  been  dropped. 

E-step:  Define  the  objective  function  to  be  the  conditional  expectation  of  the  conditional  log-likelihood  of  z  =  (y,  x),  which 
can  be  written  as 

(9) 


where  is  the  previous  estimate.  Using  ©  and  ®  and  observing  that  In  /  (x)  in  ©  is  independent  of  6,  and  hence 
irrelevant  to  the  maximization,  we  obtain 


x(Mg)  =  -Ar/2.1n(/0)--^[|y(fc)|2  +  |C'|s 

0  k=  1 

nte(y*(k)CS$(kj)] 


- 2 1 


(10) 


where  *$(*)  =  Ez{y  qU)  [3'(^)]  —  E^\y  qU)  •  Assuming  the  independence  of  csch  transmitted  symbol  snd  the 

ii 

obtain 


independence  of  x(k)  and  0]  '^ .  and  using  Bayes’  law  and  the  fact  that  ©  can  be  expressed  as  a  product  of  Ni  factors. 


X(i)(k)  —  -®x(fc)| y(k),0$  Ia;(^)] 


(0 


where 


and 


/  (x(k)  |  y(k),0$ 


c(k),d, 


/  (y(k) 


i(j) 


Pr  (x(k))  . 


f 


(k),0l$)~™  P(- 


1 


I y(k)  -  Cx(k) |2 


(ll) 


(12) 


(13) 
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M-step:  Taking  the  derivative  of  (ITOb  with  respect  to  the  real  and  imaginary  parts  of  the  complex-valued  C ,  and  then  setting 
the  results  equal  to  zero,  we  obtain  the  estimate  of  the  fading  coefficient  at  iteration  i  +  1  as 


Re  (cg}1})  =  -rr  Y, Re  (y*  x(i 

1  fc= 1 

1 

im  (£[•  j1})  =  (y*  ^ 1 


k= l 


Similarly,  maximizing  dTob  with  respect  to  the  interference  PSD  Iq  leads  to 


fU)  _ 

0,(»+l)  _  ]y 


!  M  ,  2 

1  x  ^  I  A(j)  -U), 


=  TF  E  #)-6, 


(i+l):r(i) 


fc=l 


(14) 

(15) 

(16) 


The  fading  phase  and  amplitude  can  be  explicitly  estimated  from  (IT4l)  and  (IT5]>.  but  that  is  unnecessary. 


Let  s 


U) 


=  1,2,  3, 4,  be  the  code-symbol  probabilities  obtained  from  the  soft  outputs  of  the  channel  decoder,  with 


Si  =  Pr  (x(k)  =  +1) ,  S2  =  Pr  (x(k)  =  +j) ,  S3  =  Pr  {x(k)  =  —1) ,  S4  =  Pr  (a ;(fc)  =  —j).  From  ©  and  (THU).  the  expectation 
of  x(k)  at  the  ith  EM  and  jth  receiver  iteration  is 


„(i)  Tfti)  ,  „■„(/)  n(i)  _  „(j)  TfU)  _  jjj)  p(i) 
-U)n\  61  ■“!,(*)  +  7*2  ^2,(4)  *3  ^3,(1)  7*4  -fx4,(i) 

(fc)  —  - 


"(i) 


(17) 


V  s(j) 

2=  s/3  n/3,i 

3=1 
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z[j\k)  =  log  ■ 


exp 


4J\k)  =  iog- 


exp 


exp 


JUT " 

V(*i 


—  Im  C 

*) 


/7O") 

w,-  ^2/ 

V  ‘■max  7 


JUT 

Oi(*max 
2 

r(j) 

0>(*max) 


:cW  <y* 

ax  ) 


-Rc  ( 

)  ' 

Im  ( d(tL)y *(fc)) 


exp 


exp 


j(j) 

Oi(*max) 


exp 


exp 


tO') 


— Refc^' 

/y/.  >  V  (Jmax)y 

u!(lmaxJ 


V2 


exp 


j(j) 

0i(*max) 


~JU) 

0if®max) 


■Im  c; 


<00  , 
(®max) 


(‘ 

(f tL>v'w)  +*■■ 


fU) 

0i(*max) 


hnft(iL)V’w)+vi 


(19) 


(20) 


(i) 

where  likelihood-ratio  Rp^  depends  on  the  current  CSI  estimates  as 


73  U) 
^1,(0 

=  exp 

fU)  Re(c$y(k)) 
.o,(i) 

73  GO 

=  exp 

M  C$y(k)) 

.o,(0 

3,(0 

=  exp 

fu)  Re(^fo  y(fc)) 

Jo,(0 

E>00 

=  exp 

fU) 

0,(0 


(18) 


Therefore,  for  a  given  receiver  iteration,  xjfj  (k)  and  R^\  are  updated  imax  number  of  times  using  decoder  feedback  s^.  In 
the  next  receiver  iteration,  after  channel  re-estimation,  the  fading-coefficient  and  interference-PSD  estimates  are  updated,  and 
then  used  at  the  demodulator  and  channel  decoder  to  recompute  i^+1^(fc)  and  R^J^ .  This  process  is  repeated  again  for  imax 
EM  iterations,  and  the  aforementioned  cycles  continue  likewise  for  subsequent  receiver  iterations. 

In  estimating  the  fading  parameters,  we  set  N±  =  tifb/^]  in  estimating  Iq,  we  choose  n/s  <  ufb  and  set  N\  =  71jb/2. 

£ (7)  ' 

The  EM  estimator  first  finds  the  value  of  for  a  fading  block  of  size  tifb-  Then  it  finds  the  value  of  for  each 
smaller  or  equal  interference  block  of  size  nis  using  the  value  of  found  for  the  larger  or  equal  fading  block.  When  pilot 

symbols  are  used,  we  set  xjfj  (k)  =  x(k)  for  each  known  pilot  bit,  and  there  are  no  EM  iterations  if  only  known  pilot  bits 
are  processed  in  calculating  the  channel  estimates.  The  application  of  the  EM  algorithm  is  to  obtain  both  channel-coefficient 
and  interference-PSD  estimates,  which  differs  from  m-m  where  the  emphasis  is  on  data  detection,  and  noise  statistics  are 
assumed  to  be  perfectly  known. 

Let  l  =  1,2  denote  the  two  bits  of  a  QPSK  symbol,  and  v\ .  v-j  denote  the  corresponding  log-likelihood  ratios  that  are  fed 
back  by  the  channel  decoder.  From  ED  and  l32l  Eqn.  6],  the  demodulation  metrics  (extrinsic  information)  z\^(k),l  =  1,2 
for  bits  1,  2  of  symbol  k  that  are  applied  to  the  channel  decoder  are  shown  at  the  top  of  the  next  page. 

The  number  of  EM  iterations  and  the  receiver  latency  are  reduced  by  applying  a  stopping  criterion.  Iterations  stop  once 
C^)  is  within  a  specified  fraction  of  its  value  at  the  end  of  the  previous  iteration  or  a  specified  maximum  number  is  reached. 
The  fraction  should  be  sufficiently  small  (perhaps  10%)  that  the  performance  loss  will  be  insignificant. 


IV.  BLIND  CSI  ESTIMATION 

The  EM  algorithm  in  Section  III  generates  updated  CSI  estimates  as  shown  in  (1 1 41-d  1 6b  after  the  initial  coherent  demodulation 
and  decoding  of  receiver  iteration  j  =  0.  In  [H9l-[f2Tl.  the  initial  CSI  estimates  were  obtained  with  the  aid  of  pilot  symbols. 
In  this  section,  two  methods  for  blind  estimation  of  the  initial  CSI  parameters  ^  )^o(l  ))  are  Presented, 

with  the  special  case  of  perfect  phase  information  at  the  receiver  examined  first. 


A.  Perfect  Phase  Information  at  Receiver 

The  carrier  synchronization  provided  by  a  phase-locked  loop  in  several  second  and  third-generation  cellular  standards  such  as 
IS-95  and  CDMA20001x  can  be  exploited  to  obviate  the  need  to  estimate  the  channel  phase  (which  is  also  potentially  provided 
by  2%  piloting  |27)).  Assuming  perfect  phase  information  at  the  receiver,  the  fading  amplitude  is  real-valued  and  nonnegative, 


7 


and  (fl5]>  does  not  have  to  be  computed.  A  simple  heuristic  estimate  (denoted  as  blind  method  I)  of  ),/q ^  j  for 

each  fading  block  can  be  obtained  from  the  received  symbols  as 


where 


A(0)  _  2 


nFB 


riFB  /2 

i^(fc)i 


?== i 


f  0) 

In  =  max 


D  = 


D  - 

2 


»  F  B  /  ^  o 

£„ ,  fowl2 


TlFB 


(21) 

(22) 

(23) 


represents  the  average  power  of  the  received  symbols,  and  1)  —  ,  j  is  the  difference  between  that  power  and  the 

estimated  average  power  of  a  desired  symbol.  Equation  (OTT)  would  provide  a  perfect  estimate  in  the  absence  of  noise  and 
interference.  The  parameter  h  >  0  is  chosen  such  that  j  j  1^.  .  does  not  exceed  some  maximum  value.  Ideally, 

h  is  a  function  of  Es/Nq,  but  in  this  paper  a  constant  h  =  0.1  is  always  used  for  simplicity. 


B.  Complexity  Analysis 

Although  the  EM  estimation  is  a  relatively  low-complexity  iterative  approach  to  maximum-likelihood  estimation,  it  consumes 
a  much  larger  number  of  floating-point  operations  than  pilot-assisted  schemes  do.  To  evaluate  the  complexity  of  the  EM  estimator 
in  terms  of  required  real  additions  and  multiplications  per  block  of  Ni  code  symbols,  each  complex  addition  is  equated  to  two 
real  additions,  each  complex  multiplication  is  equated  to  four  real  multiplications,  and  divisions  are  equated  with  multiplications. 
Equations  (IT4l) — (IT6l>  require  jmaximax  (6iVi  +  4)  real  additions  and  jmaxtmax  (12iVi  +  4)  real  multiplications.  Equations  (18) 
and  (19)  require  6jmaximax  real  additions,  30jmax*max  real  multiplications,  and  the  computation  of  4  exponentials.  The  initial 
estimates  calculated  using  (f2Tb — (l23l>.  which  only  need  to  be  computed  once  prior  to  the  first  EM  iterations,  require  2Ni  real 
additions,  8Ni  +  7  real  multiplications,  and  the  computation  of  the  maximum  of  two  real  numbers.  A  PACE  receiver  that  uses 
only  pilot  symbols  for  CSI  estimation  requires  6Ni  +  4  real  multiplications  and  1 2  A'|  +  4  real  multiplications  to  compute 
(fl4l>— (fl6l)  once  and  does  not  need  to  compute  the  other  equations.  Thus,  EM  estimation  increases  the  amount  of  computation 
for  CSI  estimation  by  a  factor  of  more  than  jmax?’max  relative  to  PACE. 


C.  No  Phase  Information  at  Receiver 

The  initial  CSI  estimates  proposed  in  El  and  (Ell  for  blind  method  I  are  expected  to  be  degraded  significantly  when  the 
phase  information  is  also  unknown,  since  an  arbitrary  initial  phase  value  (e.g.,  0  radians)  must  be  assumed.  To  circumvent  this 
problem,  the  initial  receiver  iteration  consists  of  hard-decision  demodulation  and  channel  decoding,  after  which  each  decoded 
bit  is  used  as  xft  Jk)  in  (HUmo.  This  step  is  followed  by  the  regular  EM  estimation  process  in  subsequent  receiver 

\  1 1 1 L  CL  X  / 

iterations.  This  approach  for  the  initial  CSI  estimates,  which  is  referred  to  as  blind  method  II  in  the  sequel,  results  in  increased 
receiver  latency  relative  to  the  previous  method  when  phase  information  is  not  available. 


D.  Blind-PACE  Estimation  Tradeoffs 

The  previously  proposed  iterative  DS-CDMA  receiver  with  PACE  lfT9ft-fi2ll  is  considered  as  the  benchmark  for  comparison 
with  the  proposed  receiver.  Assuming  an  identical  transmit-power  constraint  and  information  bit-rate  in  both  cases,  the 
elimination  of  pilots  creates  the  following  possibilities  for  methods  I  and  II: 

•  (Case  A )  An  increase  in  the  number  of  transmitted  information  symbols. 

•  (Case  B)  An  increase  in  transmitted  information-symbol  duration. 

•  (Case  C)  An  increase  in  the  number  of  transmitted  parity  symbols  (lowered  IRA  code  rate). 

The  modifications  listed  above  offset  the  loss  in  system  performance  due  to  the  degraded  CSI  estimation  obtained  from 
blind  methods  I  and  II  with  respect  to  PACE.  The  no-pilot  cases  A,  B,  and  C  have  the  same  transmitted  frame  duration  as 
the  frame  with  pilot  symbols.  Cases  A ,  B,  and  C  provide  the  most  favorable  throughput,  spectral  efficiency,  and  bit  error 
rate,  respectively.  Numerical  evaluations  of  each  of  these  cases  are  presented  in  the  next  section.  Although  a  correlated  fading 
model  is  assumed  in  the  simulations,  no  filtering  is  used  to  exploit  this  correlation  in  order  to  maintain  the  robustness  of  the 
proposed  estimator. 


V.  SIMULATION  RESULTS 


In  all  the  simulations,  the  block  sizes  are  equal,  and  the  information-bit  rate  is  100  kb/s.  Increasing  the  block  sizes  increases 
the  accuracy  of  the  EM  estimators,  but  decreasing  the  block  sizes  allows  closer  tracking  of  the  channel  parameters  and 
includes  more  diversity  in  the  receiver  computations.  In  most  of  the  simulations,  except  where  stated,  we  set  nFi=npB  =  40 
and  spreading  factor  g  =  31.  The  number  of  closed-loop  receiver  iterations  is  set  to  jmax  =  9,  as  there  is  insignificant 
performance  improvement  for  jmax  >  9.  The  number  of  internal  EM  iterations  is  zrnax  =  10.  There  is  one  decoder  iteration 
per  receiver  iteration.  A  IRA  code  (data  block  size  K  =  1000)  with  sum-product  algorithm  decoding  fT8ll  is  used  without 
channel  interleaving.  The  IRA  code  is  rate-1/2  when  PACE  is  used.  Jakes  correlated  fading  of  the  desired  signal  and  a  mobile 
velocity  of  120  km/hr  are  assumed.  Flat  fading  is  assumed  in  most  of  the  simulations,  whereas  a  frequency-selective  channel 
is  examined  in  Section  IV-EI  The  iterative  PACE  receiver  considered  for  comparison  contains  9.1%  pilot-symbol  overhead, 
which  has  been  shown  to  have  a  decoding  performance  close  to  the  conventional  3GPP  LTE  receiver  EQ.  For  each  scenario 
tested,  5000  Monte  Carlo  simulation  trials  were  conducted.  To  avoid  repetition,  a  selection  of  representative  examples  out 
of  the  many  possible  combinations  of  channel  coding,  phase  information,  interference  models,  and  no-pilot  modifications  are 
presented  next. 

The  bit  error  rate  (BER)  is  calculated  as  a  function  of  E^/Nq,  where  Ei,  =  (N/2K)ES  is  the  energy  per  bit.  The  information 
throughput  is  a  vital  performance  criterion  in  addition  to  the  BER.  One  of  the  primary  motivations  in  removing  pilot  symbols 
is  the  expectation  of  achieving  greater  throughput,  even  though  the  BER  performance  may  be  degraded  marginally.  We  define 
throughput  R  as 


R  = 


information  bits  in  a  codeword 
codeword  duration 


x  (1  —  BER)  bits/s. 


(24) 


A.  Single-user  environment,  perfect  phase  knowledge 


Fig.  3.  BER  versus  E^/Nq  for  IRA-coded  iterative  receiver  in  single-user  environment  with  phase  provided  by  PLL. 
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Fig.  4.  Information  throughput  versus  E^/Nq  for  IRA-coded  iterative  receiver  in  single-user  environment  with  phase  provided  by  PLL. 


For  the  first  set  of  results  in  Figs.  00  a  single-user  environment  and  perfect  phase  knowledge  at  the  receiver  are  assumed. 
Fig.  0  displays  the  BER  versus  E^/Nq  for  an  IRA-coded  iterative  receiver  operating  with  perfect  CSI,  PACE,  blind  method  I 
with  cases  A,  B,  and  C ,  and  blind  method  II  with  cases  A  and  C,  respectively.  The  key  observation  is  that  blind  method  II 
is  worse  than  method  I  by  2  dB  at  BER  =  10-3  for  both  case  A  and  case  C,  which  illustrates  the  well-known  sensitivity  of 
the  EM  algorithm  to  the  accuracy  of  the  initial  estimates. 

The  addition  of  extra  parity  bits  to  blind  method  I  (case  C,  rate- 1000/2200)  offers  the  greatest  improvement  in  BER, 
surpassing  even  the  rate-1/2  code  with  perfect  CSI  at  high  E^/Nq.  The  increase  in  number  of  information  symbols  (case 
A)  results  in  the  worst  BER  performance  with  a  separation  of  1  dB  and  0.5  dB  from  PACE  and  case  B  at  BER  =  10-3, 
respectively. 

The  various  scenarios  featured  in  Fig.  0  were  also  tested  under  a  slow-fading  channel  with  mobile  velocity  of  10  km/hr.  It 
was  observed  that  all  the  BER  curves  were  shifted  towards  the  right  by  up  to  7  dB  at  BER  =  10“3,  but  the  overall  trends 
among  the  different  cases  remained  the  same. 

Fig- IH exhibits  information  throughput  R  versus  A/, / Ar(l  for  the  IRA-coded  iterative  receiver  with  the  scenarios  of  Fig.0  The 
throughput  advantage  of  case  A  is  achieved  even  though  no  pilot  symbols  are  used  at  all;  i.e.,  the  initial  estimation  is  totally 
blind.  It  is  evident  that  increasing  the  symbol  duration  or  adding  additional  parity  information  does  not  give  the  proposed  blind 
methods  any  significant  advantage  in  throughput  over  PACE.  Both  blind  methods  with  cases  B,C  and  PACE  provide  about 
20%  less  throughput  than  the  receiver  with  perfect  CSI. 

B.  Multiuser  environment,  unknown  phase 

A  4-user  interference  environment  with  equal  mean  bit  energies  for  all  users  at  the  receiver,  E^/Nq  =  20  dB,  and  no 
phase  information  at  the  receiver  is  examined  next.  It  is  assumed  that  both  the  interference  levels  and  the  unknown  phase  are 
constant  during  each  subframe.  Each  interference  signal  experiences  independent  Jakes  correlated  fading  and  uses  independent 
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data  and  Gold  sequences  with  respect  to  the  desired  signal.  The  simulation  uses  chip-synchronous  interference  signals,  which 
is  a  pessimistic  worst-case  assumption  El.  Two  variations  of  CSI  estimation  are  examined  here:  partially  adaptive  with 
only  fading  coefficient  C^)  being  estimated  using  (ITdl).  (fl5l).  and  set  equal  to  No  for  all  subframes;  and  fully  adaptive 

estimation  of  both  and  using  (fl4t.  (IT5l).  and  (fl6l>. 
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Fig.  5.  BER  versus  E^/Nq  for  IRA-coded  iterative  receiver  affected  by  MAI  from  4  users,  fully  and  partially  adaptive  estimation,  and  unknown  phase. 

Fig.  [5] displays  IRA-coded  BER  versus  If,/ No  for  partially  and  fully  adaptive  CSI  estimation  per  fading  block  and  case  C 
for  both  blind  methods.  The  mismatch  of  Iq  and  the  true  value  of  /(l  at  the  demodulator  and  decoder  results  in  a  high  error 
floor  for  the  partially  adaptive  cases.  The  intuition  behind  the  error  floor  is  that  the  partially  adaptive  estimator  overestimates 
the  true  signal-to-interference-plus-noise  ratio  (SINR)  by  disregarding  the  MAI,  with  the  degree  of  overestimation  increasing 
with  SINR.  For  IRA  codes,  it  was  shown  in  l33ll  that  both  under-  and  overestimation  of  the  SINR  degrades  the  IRA  decoder 
performance.  The  fully  adaptive  estimation  offers  a  more  accurate  SINR  estimate  and,  hence,  suppresses  interference  and 
reduces  the  error  floor  significantly.  This  interference  suppression  is  achieved  without  using  the  far  more  elaborate  multiuser 
and  signal  cancellation  methods  that  could  be  implemented  in  a  DS-CDMA  receiver.  For  both  partially  and  fully  adaptive 
estimation,  it  is  observed  that  blind  method  II  now  outperforms  method  I  due  to  better  phase  estimation,  whereas  both  blind 
methods  outperform  PACE  at  BER  =  10-3  due  to  the  added  parity  information. 

Fig- E] demonstrates  the  IRA-coded  receiver  throughput  offered  by  the  proposed  methods  under  MAI  from  4  users.  The  blind 
methods  always  provide  a  better  throughput  compared  with  PACE;  for  example,  method  I  with  case  A  is  superior  by  9%  to 
both  PACE  scenarios  when  Eb/No  >  5  dB.  It  is  observed  that  both  partial  and  fully-adaptive  estimation  methods  offer  a  similar 
asymptotic  throughput,  which  indicates  that  partial  CSI  estimation  may  be  sufficient  for  applications  with  a  non-stringent  BER 
criterion.  On  the  other  hand,  error-critical  applications  requiring  less  than  BER  =  1  (V  3  must  use  the  fully  adaptive  CSI 
estimation,  as  seen  from  Fig.  [5] 
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Fig.  6.  Information  throughput  versus  E^/Nq  for  IRA-coded  iterative  receiver  affected  by  MAI  from  4  users,  fully  and  partially  adaptive  estimation,  and 
unknown  phase. 


C.  Varying  fading-block  size ,  unknown  phase 

In  urban  mobile  environments,  the  phase  can  be  expected  to  change  significantly  after  approximately  s  to  '-jA  s,  where 
fd  is  the  maximum  Doppler  shift.  For  the  assumed  mobile  velocity  of  120  km/hr,  this  time  range  corresponds  to  roughly  10  to 
40  code  bits  at  100  kb/s.  The  fading  and  interference  block  sizes  ufb  —  niB  are  therefore  varied  accordingly,  and  no  phase 
information  is  assumed  to  be  available  at  the  receiver  for  the  next  set  of  results. 

Fig.  [7]  displays  fully  adaptive  IRA-coded  BER  versus  E^/Nq  for  blind  methods  I  and  II  with  case  G,  9.1  %  PACE,  and 
perfect  CSI  decoding  for  ufb  =  10  and  40  in  a  single-user  environment.  An  improvement  of  1  to  2  dB  was  observed  for  all 
methods  for  the  smaller  fading-block  size  of  ufb  =  10  due  to  the  increased  fading  diversity.  The  throughput  with  case  A  is 
shown  in  Fig.  [8]  It  is  observed  that  the  throughput  gains  of  the  proposed  blind  methods  over  PACE  (roughly  9%  at  intermediate 
to  high  Eb/No)  are  preserved  even  when  the  phase  is  initially  unknown  at  the  receiver. 

D.  Varying  MAI,  unknown  phase 

IRA-coded  iterative  receiver  performance  with  blind  method  II,  case  C  is  examined  for  3  and  6  MAI  signals  with  equal  mean 
bit  energies  for  all  users  at  the  receiver  in  Fig.  [9]  The  partially  adaptive  estimation  is  unable  to  cope  with  the  interference  caused 
by  6  MAI  signals  regardless  of  the  spreading  factor,  whereas  the  fully  adaptive  estimation  offers  a  substantial  improvement  in 
BER.  The  benefit  of  an  increased  spreading  factor  (g  =  127  versus  g  =  31)  is  more  apparent  at  low  bit  error  rates  for  fully 
adaptive  estimation.  For  example,  the  fully  adaptive  estimation  with  3  MAI  signals  improves  by  a  factor  of  approximately  5 
dB  at  BER  =  10-5,  despite  nonorthogonal  spreading  sequences  and  imperfect  CSI. 
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Fig.  7.  BER  versus  E^/Nq  for  IRA-coded  iterative  receiver  in  single-user  environment,  varying  npg,  and  unknown  phase. 


E.  Multipath  channel 

A  DS-CDMA  system  can  exploit  a  frequency-selective  fading  channel  by  using  a  Rake  receiver.  As  an  example,  we  assume  a 
channel  with  three  resolvable  multipath  components  (with  known  delays)  of  the  desired  signal  and  a  Rake  combiner  with  three 
corresponding  fingers.  The  multipath  components  undergo  independent  fading  across  the  fingers,  but  follow  the  Jakes  correlated 
fading  assumption  over  time.  The  multipath  components  follow  an  exponentially  decaying  power  profile  across  the  fingers, 
i.e.,  E  [ ai ]2  =  e_(z_1),  1  =  1,2, 3.  Each  interference  signal  has  the  same  power  level  in  each  finger  and  undergoes  independent 
Jakes  correlated  fading.  The  assumption  of  independent  multipath  fading  amplitude  and  phase  coefficients  for  the  desired  signal 
allows  us  to  apply  the  proposed  EM-based  channel  estimation  scheme  separately  in  each  finger.  The  Rake  combiner  performs 
maximal-ratio  combining  (MRC)  of  the  received  symbol  copies  based  on  channel  and  interference-PSD  estimates  computed 
at  all  fingers.  The  MRC  decision  statistic  obtained  from  the  Rake  combiner  is  then  passed  to  the  QPSK  demodulator  metric 
generator,  which  generates  soft  inputs  for  the  common  channel  decoder.  The  channel  decoder  soft  outputs  are  fed  back  to  the 
three  channel  estimator  blocks,  which  then  recompute  updated  channel  coefficients,  as  described  in  Section  III. 

Fig.  [TO]  displays  the  Rake  receiver  performance  for  various  levels  of  MAI  with  Method  II  under  case  C,  where  all  users 
have  length-127  Gold  sequences.  It  is  observed  that  the  additional  diversity  due  to  Rake  combining  improves  performance  as 
expected,  but  the  performance  disparity  between  partially  and  fully  adaptive  estimation  remains  large. 

VI.  Conclusions 

It  has  been  shown  that  pilot  symbols  are  not  essential  to  the  effectiveness  of  DS-CDMA  receivers  with  coding,  coherent 
detection,  and  channel  estimation.  If  the  pilot  symbols  are  replaced  by  information  symbols,  the  throughput  increases  relative 
to  PACE  whether  or  not  interference  is  present.  If  the  BER  is  the  primary  performance  criterion,  then  replacing  the  pilot 
symbols  by  parity  symbols  gives  a  lower  BER  than  PACE.  If  the  spectral  efficiency  is  of  primary  importance,  then  extending 
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Fig.  8.  Information  throughput  versus  E},/Nq  for  IRA-coded  iterative  receiver  in  single-user  environment,  varying  npg,  and  unknown  phase. 


the  symbol  duration  after  the  removal  of  the  pilot  symbols  offers  an  improvement  relative  to  PACE,  albeit  at  the  cost  of  a 
slight  increase  in  the  BER. 

The  estimation  of  the  interference  PSD  has  been  shown  to  enable  the  significant  suppression  of  interference.  This  suppression 
is  achieved  without  using  the  far  more  elaborate  multiuser  and  signal  cancellation  methods  that  could  be  implemented  in  a 
DS-CDMA  receiver. 
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